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Self- organized monolayers of highly flexible Frechet dendrons were deposited on graphite surfaces 
by solution casting. Scanning tunneling microscopy (STM) reveals an unprecedented variety of 
patterns with up to seven stable hierarchical ordering motifs serving as a versatile model system. The 
essential molecular properties determined by molecular mechanics simulations are condensed to a 
coarse grained interaction site model of various chain configurations. In a Monte Carlo approach with 
random starting configurations the experimental pattern diversity can be reproduced in all facets 
of the local and global ordering. Based on an energy analysis of the Monte Carlo and molecular 
mechanics modeling the thermodynamically most stable pattern is predicted coinciding with the 
pattern, which dominates in the STM images after several hours or upon moderate heating. 

PACS numbers: 68.43.-h, 81.16.-c, 68.37.Ef 



Organic molecules mainly self-organize via hydro- 
gen bonds, [l|, m metal-coordination |3i] and van der 
Waals interactions Voices have been raised, [HI, [6| 
that more innovative modeling approaches 0, Hi are 
necessary in order to analyze and predict molecular 
patterns for fostering applications in surface functional- 
ization sensors (loll , catalysis [11], and in molecular 
electronic devices [12|. Predictability could speed up 
molecular design and advance the understanding of the 
self-organization process itself. For predicting several 
coexisting ordering motifs (cf. 13|, a new multi- 



modeling approach is required, that does not r ely on the 
prior knowledge of the resulting patterns |15Ml8l|. Here 
we show, that based on a coarse grained interaction site 
model (omitting chemical details) we are able to not only 
reproduce and predict all features of the local/global 
ordering motifs of self-organized molecular layers, but 
also independently predict the thermodynamically most 
stable pattern and thus deliver an innovative approach 
to the understanding of molecular self-organization. 

In this article, the self-organization of flexible supra- 
molecular building blocks, Frechet dendrons [H, |2q| . 
is investigated in highly resolved scanning tunneling 
micrographs and modeled with various techniques: 
density functional theory (DFT), molecular mechanics 
(MM) and Monte Carlo (MC) simulation based on 
an interaction site model. The second generation 
asymmetric Frechet dendrons consist of three phenyl 
rings, two of which are elongated by two pairs of alkoxy 
chains, a longer 12mer pair and a shorter 8mer pair (see 
Fig. [1] left), responsible for van der Waals and weak 



hydrogen bonds (28|. All van der Waals interactions 
present are the sum of the intermolecular interaction 
(chain - neighboring - chain) and the molecule-substrate 
interaction (molecule - substrate). By varying the 
internal angles of the pairs of alkoxy chains in an 
interaction site model, all facets of the experimentally 
found patterns can be predicted within our approach. A 
zero temperature analysis of the MC simulations points 
to a thermodynamically most stable pattern in line with 
conclusions of the local energetics derived from MM 
simulations and experimental findings upon heating. 
For the STM measurements a Nanoscope Multimode III 
equipped with a low-current converter under ambient 
conditions is employed. Monolayers of Frechet dendrons 
( 3 , 5-bis [ (3 , 5-bisocty loxypheny 1) met hy loxy] phenyl esther , 
synthesis see [21]), are cast from 0.2 mM solutions 
in ethanol, hexane or heptadecane on highly oriented 
pyrolytic graphite (HOPG) surfaces. The emerging pat- 
terns are analyzed by combining several complementary 
approaches: a.) DFT, utilizing the CASTEP module 
of Material Studio 4.3 employing Perdew-Wang '91 
generalized grad ient approximation exchange correlation 
functionals (22| and a plane wave basis set with an 
energy cutoff at 260 eV; b.) MM, performed on a fixed 
double-layer of graphite (universal force fields (23[ in 
periodic boundary conditions) based on starting values 
obtained from the STM experiments; c.) MC annealing 
simulations in Metropolis scheme based on 60° rotations 
of a coarse grained interaction site model, using periodic 
boundary conditions for minimizing finite-size effects 
in system sizes ranging from several hundreds to a few 




FIG. 1: (color online) The dode- 
cyl/octyl terminated Frechet dendron 
is highlighted in yellow. One of seven 
different ordering patterns on HOPG 
is exemplified: middle: STM image of 
a jigsaw pattern of the Frechet den- 
dron (12.5 nm x 7 nm, Usias = - 
800 mV, |/t| = 8 pA). left: Atomistic 
modeling by an energy minimized MM 
simulation (see text), right: simula- 
tion of the LDOS by DFT calculations 
based on the input of the geometry de- 
termined by MM. 



FIG. 2: (color online) (a-d) STM results and corresponding (e-f) molecu- 
lar mechanics modeling (see text) of four experimentally found patterns of 
Frechet dendrons : a) sawtooth b) trimer honeycomb, c) jigsaw, d) tiretrack. 
The molecular backbones marked in yellow indicate the local ordering motif. 
(12.5 nm x 22 nm, Usias = -700 to -800 mV, |/t| = 5 to 30 pA) 



thousand objects. 

The mixed graph of Fig. [1] ^ exemplifies MM and DFT 
simulations together with measurements of the jigsaw 
pattern. After applying a droplet of the molecules in 
hexane, the jigsaw pattern {p2 symmetry) appears 
as one of seven general ordering motifs in about 25% of 
the substrate area. A detailed analysis of many STM 
images yields starting configurations jsoj for a MM 
energy minimization [31,] to calculate molecular confor- 
mations (see left of Fig. [T]) and relative positions. From 
the MM determined adsorption energy the dominant 
van der Waals part ^] can be extracted, which will be 
discussed later. 

As the energy-minimized MM simulations can only 
provide atomic positions and not electronic states, we 
used these energy minimized geometries as input to 
derive the local density of states (LDOS) of a free single 
molecule usine; DFT. Convoluted density contours of 
a planar slice in the LDOS are plotted fs^ (right, Fig.[T]) 
for a direct comparison with the STM measurement: 
the phenyl rings of the molecular cores are discernible as 
three bright protrusions and atoms of the alkoxy arms 
emerge (visible in the dark areas) as series of faint lines 
(right. Fig. [1]) in agreement with the experimental STM 
images (middle. Fig. [1]). 

In a large STM study, we have imaged seven ordering 
motifs of this Frechet dendron in excellent resolution. 
Here we focus on four main ordering motifs: sawtooth 
: trimer honeycomb : jigsaw : tiretrack (see Fig. [2|). 
20 nm x 20 nm sized domains of these patterns in ratios 
1% : 15% : 25% : 20% typically cover an HOPG surface 
half an hour after applying a hexane droplet. About 39% 
of the HOPG surface is occupied by domain boundaries, 
mobile molecules or two other patterns: upright and 



small tiretrack. Minutes after casting a droplet of the 
Frechet dendron in hexane, the honeycomb pattern 
dominates. The sawtooth pattern seldom appears as an 
independent pattern and often as a stacking fault of the 
honeycomb pattern. After several hours typically 60% of 
the HOPG surface is covered with the tiretrack pattern 
in hexane. 

As DFT calculations can only poorly describe van der 
Waals interactions, the dominant part of the internal 
energy, the adsorption energy per unit cell or area, is 
evaluated from the MM minimization energy f35|. This 
covers the dominant intermolecular interactions stem- 
ming from the alkoxy chains, statistically contributing 
an interaction strength of 14 kJ/mol per four CH2-units 
for more than half but not fully interdigitated alkoxy 
chains on graphite (previous studies yield 18 kJ/mol 
for fully and 8.5 kJ/mol for half-interdigitated alkoxy 
chains on graphite surfaces (lH). We have calculated 
the van der Waals part of the total energy I) of a 
substrate supported monolayer in MM energy minimized 
geometry, isup^ H) of a gas phase net [2], ig^ and 
HI) of one isolated molecule, iiso- With that we can 
carefully separate the intermolecular interaction (chain - 
neighboring-chain) part /c-n-c/molecule= (ign/n) — iiso 
from the substrate interaction (molecule - substrate) 
part _ 5 /molecule = iign — Uup)/^ of the van der Waals 
energy per molecule, with n, the number of molecules 
per unit cell (see Table I a). With dividing accordingly 
by the area per unit cell, the respective energies per nm^ 
are derived (see Table I b). 

Molecular domains grow or dissolve by outer molecules, 
which have only a fraction ^ < 1 of the total number 
of neighboring molecules, adsorbing or going over into 
a mobile state, respectively. When initially forming 
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TABLE I: Van der Waals part of the adsorption energies de- 
termined by MM energy minimization, a) per molecule, b) 
per area. The labels refer to: m-s, molecule-substrate inter- 
actions and c-n-c, chain-neighboring-chain interactions. 
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FIG. 3: (color online) Top: angle configuration for conformer 
a. The box displays the six possible orientations on a fixed 
lattice of hexagonal symmetry. Bottom: Conformers /3, 7, 
and 5 with the respective internal angles indicated. 



domains after applying a droplet, the gain in van der 
Waals adsorption energy Im-s + ^{Ic-n-c) per molecule 
rules the process: the honeycomb pattern in table I a. 
Over time the HOPG surface is covered with some 
100 nm^ sized domains, which transform according to 
the Kitaigorodskii principle of avoiding free space [26|. 
Hence the highest gain in van der Waals adsorption 
energy Im-s + I c-n-c P^'^ ctrea dominates the outcome: 
the tiretrack pattern in table I b. This agrees with 
the thermodynamically stable pattern experimentally 
identified by moderate heating, which is the tiretrack 
pattern in hexane. Thus, the initial domains form open 
pore patterns dominated by the highest energy gain per 
molecule, while with time the pattern with the largest 
energy gain per area wins. However, all calculated 
energies are very close, reflecting that the phases can 
coexist. 

In order to identify the ground state energy we 
performed independent MC simulations, based on coarse 
graining the dominant molecule and substrate properties 
in an interaction site model. Rather than including as 
much chemical detail as possible, we investigated the 
crucial properties reproducing key experimental findings. 
When large molecules adsorb on ten or more substrate 
atoms the organization is dominated by steric hindrance; 
the role of the substrate is to align the molecular 
backbone and to some extend the molecular chains. In 
this work, we systematically varied the internal angles 
of the pairs of longer and shorter arms of interaction 
centers (see Fig. [3]), thereby probing different molecular 
conformations, in order to determine how critical the 
angles are on the predictive power of the interaction site 
model. 



As in our previous work flFl the substrate symmetry 
is accounted for by allowing discrete 7r/3 rotations of 
the coarse grained molecules, see Fig. [3] a, on fixed 
sites in a lattice of hexagonal symmetry. The geometry 
and spacing of the molecular backbone is transferred 
from the MM modeling. The chains are modeled by 
a small number of beads, disposed in rigid, straight 
arms at distances determined by MM and in |15| 
approximate angles of 100° and 145° between long and 
short pairs of alkoxy chains, respectively. We refer to 
a Frechet -dendron with the latter chain conformation 
as conformer a in this article. Here, the internal angle 
within each pair of arms has been systematically varied 
from 180° and 60°, referred to as conformer /3, over 120° 
and 120°, conformer 7, to 60° and 180°, conformer 5, for 
long and short pairs of alkoxy chains, respectively. One 
bead represents four CH2-units, exhibiting short-range 
van-der- Waals attractions, which are described by a 
Lennard- Jones potential, U{r) = 4:e[{a/r)^'^ — (cr/r)^]. 

The MC simulations were performed by preparing 
random starting configurations (see o) in Fig. 2]) and 
relaxing the system by slowly "cooling". The lattice 
constant a, accounting for the density, has been varied 
between 2.95 a and 4.0 <j, corresponding to a change in 
coverage density by a factor of 2.5, while the dimen- 
sionless inverse temperature e/ksT ranged from 0.025 
to 6. The ordered patterns obtained: i), ii), iii), iv) 
in Fig. [4j correspond well in long-range symmetry and 
local ordering with the structures found experimentally 
a), b), c), d) in Fig. [21 respectively. While the chain 
configuration of conformer a (100° and 145° for long 
and short pairs of alkoxy chains, respectively) leads to 
a striking variety of patterns (see [15|), the systematic 
choice of internal angles in conformers /3, S and 7 
narrows the pattern variety. 
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ordering states. Though, ah three honeycomb patterns 
have a similar appearance, upon closer examination, 
only the trimer honeycomb pattern ii) displayed in 
Fig. m correctly describes the experimentally found 
pattern b) of Fig. [2l Hence, the honeycomb and inverted 
honeycomb pattern find no explicit resemblance in the 
experimental data and are unstable in a MC simulation 
of finite temperature. 

In contrast to patterns generated from random starting 




FIG. 4: (color online) Monte Carlo (MC) simulations (i- 
iv). 0) Random starting configurations allowing discrete n/3 
rotations of the coarse grained molecules. Upon slow cooling 
the following patterns formed from random state (0) indicated 
with red arrows: i) Saw tooth pattern for conformer /? and 
for conformer a, ii) trimer honeycomb pattern for conformer 
(5, iii) jigsaw pattern for conformer a, iv) tiretrack pattern for 
conformers a, j3, and 7. All phases were stable upon heating 
to 300° C. 
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FIG. 5: (color online) Alternative phases predicted condensed 
from random starting configurations by the Monte Carlo ap- 
proach:, but not found in the experiment: v) honeycomb pat- 
tern of conformer a and 7, and vi) inverted honeycomb pat- 
tern of conformer 7, can be condensed from the random start- 
ing configuration. Both the inverted honeycomb and honey- 
comb MC patterns proved unstable upon simulating a finite 
temperature and thus cannot be considered as likely patterns; 
these patterns find no close resemblance in the experimental 
data. 



The patterns that can be generated from the respective 
conformers are: sawtooth and tiretrack for j3 (see i) 
and iv) in Fig. Hj), tiretrack, honeycomb and inverted 
honeycomb for 5 (see iv) in Fig. [H and v) and vi) in 
Fig. [5l respectively) and trimer honeycomb for 7 (see 
ii) in Fig. Hj). So with the latter conformer 7 even the 
intriguing open-pore host structure of molecular trimers 
arranged in hexagonal symmetry (see b) in Fig. [2]) can 
now be generated from random starting configurations 
by slow cooling. The trimer honeycomb pattern b) in 
Fig. [2] displays an organizational chirality with an in- 
cluded guest molecule of non-fixed orientation (described 
in another experiment (l^l). Hence all experimentally 
found patterns can be completely described. The MC 
patterns i) to iv), in very good correspondence with 
experimental findings, proved to be reliable, as raising 
the temperature showed that these phases stay ordered 
up to 300° C and can hence be considered as stable 
states. 

Additionally, we exemplify here two more patterns (see 
Fig. [5|) named v) honeycomb pattern and vi) inverted 
honeycomb pattern (both found for conformer 7), that 
were condensed from random starting configurations, 
but could not be verified in the STM experiments. The 
stability of these MC patterns was tested by simulating 
a finite temperature; both patterns proved unstable 
upon "heating" and thus cannot be considered likely 



configurations, patterns i), ii), iii), iv) (see Fig. [3]) can 
be pre-prepared for all conformers allowing to compute 
the associated energies in a zero temperature energy 
analysis. So the state of lowest energy at T = can 
be identified in each case as function of the density. 
In the case of artificially prepared patterns the energy 
analysis can only indicate how energetically attractive 
a structure is, not if the pattern is likely. The pattern 
that minimizes the energy varies, depending on the 
packing fraction (see bar graphs in Fig. [6] a) and b) 
and the internal angle, hence, the conformer. In the 
experimentally relevant regime (3.2a/cr), the most 
favorable configuration for conformer a (see [lH) and 
conformer j3) is the tiretrack pattern (see Fig. [6] a). 
As this pre-prepared pattern is stable upon "heating" 
to 300° C and in a second way can be generated from 
random starting configurations, the tiretrack pattern is 
identified as the thermodynamically most stable pattern 
in nice agreement with the experimental findings and 
the MM energy analysis discussed above. 
For the chain configurations of conformer 5 and 7 the 
sawtooth and trimer honeycomb patterns (Fig. [6] b)) 
appear to have the lowest energy at T = in the 
experimentally relevant regime, respectively. Probing 
the stability of the latter two patterns by raising the 
temperature leads to a destabilization and thus clearly 
indicates that neither the sawtooth nor trimer honey- 
comb pattern can represent the thermodynamically most 
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stable pattern around room temperature. 
In summary, the systematic variation of internal angles 
of the pairs of chains (representing the chain-chain 
interdigitation) reveals that the local and global order- 
ing motif of all experimentally found patterns can be 
reproduced by cooling random starting configurations. 
The zero temperature analysis of pre-prepared patterns 
in all combinations of conformers/patterns underlines 
that only the tiretrack pattern found for conformer a 
(see [lH) and (5 represents an ordered state with lowest 
energy at T = and is at the same time stable upon 
"heating". Thus, the tiretrack pattern denotes experi- 
mentally and theoretically the thermodynamically most 
stable motif at finite temperatures, even for different 
molecular conformations in the interaction site model. 
Yet, patterns cannot only be reproduced by the MC 
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FIG. 6: (color online) conformer /3, coarse grained model 
with angles 180° and 60° between the long and short chains; 
a.) zero temperature energy analysis for determining the as- 
sociated lowest energy state; inset: triangular lattice display- 
ing the six possible orientations . The experimental relevant 
regime corresponds to 3.2a/cr. b.) The pattern of lowest 
energy in the zero temperature analysis for conformations 7 
and 5. At the experimentally relevant regime of 3.2a/cr, the 
sawtooth pattern for conformer 7 and jigsaw pattern for con- 
former (5, did not prove stable upon heating and hence cannot 
represent the thermodynamically most stable pattern. 



(Fig.©. 

Our study underlines that main geometrical features 
are the dominant driving force of molecular self- 
organization. The various patterns found experimentally 
for the Frechet dendron presented in this study serve 
as a versatile model system to systematically test our 
multi-modeling approach. Five structural phases have 
been closely modeled by MM energy minimization; their 
quality is underlined exemplarily by a DFT LDOS sim- 
ulation plotted adjacent to the STM image. Evaluating 
the van der Waals part of the MM determined adsorption 
energy favors an open pore trimer honeycomb structure 
in the initial phase of adsorption (energy per molecule 
maximized) and a densely packed pattern reminiscent 
of tire tracks after various phase transformations in the 
final state (energy per area maximized). The calculated 
energies are very close reflecting the coexistence of 
phases for competing ordering motifs. Moreover, we 
have refined our independent MC approach with a coarse 
grained interaction site model of various internal angles 
and so successfully reproduce the experimental pattern 
diversity in all facets of the local and global ordering. 
The lowest energy pattern at T = can be identified in a 
zero temperature analysis for the different chain confor- 
mations. All theoretical and experimental findings point 
to the tiretrack pattern as thermodynamically most 
stable pattern for the dodecyl/octyl terminated Frechet 
dendron. Furthermore, a newly predicted local ordering 
motif has been experimentally verified afterwards. We 
thoroughly and successfully tested an interaction site 
approach to molecular self-organization, which will serve 
as an innovative analysis tool in the future fostering the 
application of self-organized molecular monolayers in 
various areas of science. 
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FIG. 7: (color online) a.) Another local ordering motif of 
the MC simulation Fig.|4]iv.) found afterwards by b.) STM: 
measurement of a rarely found fifth pattern with correspon- 
ding b.) MM modeling.( 12.5 nm x 22 nm, UBias = -800 mV, 
\It\ = 8 pA) 



approach, but also predicted. In Fig. S] iv) a second 
local ordering motif with a side by side arrangement of 
the molecular backbones is apparent (Fig. [7] a). This 
local ordering motif named "wave" has been discovered 
after the theoretical prediction; it occurs rarely for 
dodecyl/octyl terminated Frechet dendrons in ethanol 
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